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Abstract 

We study higher-order effects in the electromagnetic production of electron- 
positron pairs in relativistic heavy ion collisions. Treating the field of the 
heavy ions as an external field and neglecting the interaction among electrons 
and positrons, we show that the A^-pair creation amplitude is the antisym- 
metrised product of one-pair creation amplitudes and the vacuum ampli- 
tude. Neglecting contributions coming from exchange terms, we show that 
the total probability for N pairs is approximately a Poisson distribution. We 
investigate further the structure of the reduced one-pair amplitude, concen- 
trating especially on multiple-particle corrections. We calculate the first of 
these corrections in second order Magnus theory based on our previous result 
in second-order Born approximation for impact parameter b zero. Explicit 
calculations show that the total probability is increased up to 10 % by this 
correction for realistic collider parameters. The calculations can also be used 
to confirm the use of the Poisson distribution for the total probability. 

12.20.-m, 34.90.-hq, 14.60.Cd 
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I. INTRODUCTION 



The electromagnetic production of electron-positron pairs in relativistic heavy ion col- 
lisions has gained some interest recently due to the observation that the total probability 
violates unitarity if calculated in lowest order even for realistic energies and impact parame- 
ters as large as the Compton wavelength [|l|. It can be expected therefore that higher-order 
effects — especially the multiple-pair production — are of importance. Several models have 
been used in order to cure the unitarity violation and predict multiple-pair probabilities, 
based on the summation of some subclasses of diagrams All authors find that the 

pair multiplicity can be approximated by a Poisson distribution and that the lowest-order 
results should be interpreted as the pair multiplicity. Unitarity is no longer violated due to 
the inclusion of the vacuum amplitude, which reduces all A^-pair probabilities by the same 
factor. Recently lonescu has proposed a model, which differs from the previous ones, but 
which is not applicable for multiple-pair production. 

We are going to extend these existing models into two directions: First we address 
the question how the A^-pair creation amplitude is exactly reducible to one-pair creation 
amplitudes and what approximations are necessary in order to get the Poisson distribution. 
We then compare with earlier results. Secondly we study higher-order processes contributing 
to the lowest-order one-pair creation amplitude. We concentrate on multiple-particle effects, 
which have not been included in previous models. 

This paper is arranged as follows: In Sec. |I| we show that the A^-pair creation amplitude 
can be written as an antisymmetrised product of reduced one-pair creation amplitudes and 
the vacuum amplitude. For this we use Feynman boundary conditions for the fermion field. 
Neglecting contributions from exchange terms in the calculation of the total probability 
we get a Poisson distribution for the total A^-pair probability quite easily. In Sec. |T| the 
same result is derived using the perturbation expansion of the ^-matrix, which gives us 
some insight into the processes contributing to the reduced pair creation amplitude and the 
vacuum amplitude. Finally we compare in Sec. |V| our results with other existing ones. 

In Sec. 0we derive the general form of the S'-operator in second-order Magnus theory and 
use it in Sec. [V| in order to extract from this the lowest-order multiple-particle correction to 
the one-pair creation probability. This correction is then calculated in Sec. |VI1| based on our 
analytic form of the differential pair creation amplitude in second-order Born approximation 
for impact parameter b zero [0 . The same calculation is used also to study the deviation of 
the total two-pair creation probability from the Poisson distribution, which is found to be 
rather small. Results of these calculations together with conclusions are then summarized 
in Sec. [Vmt 

In the Appendix, we show a way to calculate the coefficients appearing in the Magnus 
theory expansion, which can therefore be calculated to arbitrary order. 

Throughout this paper, we treat the fields of the heavy ions as external fields. We neglect 
the electromagnetic interaction between electrons and positrons. This seems to be justified 
for almost all fermions, apart from electrons and positrons, which are produced with very 
low relative velocitiy, and even in a bound e~^-e~ state, due to the smallness of the coupling 
constant a between the fermions compared to the effective coupling constant Za for the 
interaction with the external field. Due to these assumptions, our many particle system is 
essentially a system of independent fermions in an external field. For such a system the S 
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operator is known to be of the form of a time-ordered exponential 



S = Texp|-i J Hj{t)dtj 

= r exp jy" d^x : ^{x)[-ie 4{x)]'^{x) :| , (1) 

where Hj is the Hamiltonian in the interaction picture. This form of the S operator is the 
starting point for our calculations. 

The general theory of a fermion interacting with an external electromagnetic field is 
rather old and has been developed already from the beginning of QED. It has been developed 
especially by Feynman |Q and Schwinger 0. An overview can also be found in [p!0|-[T2|. 

Although we discuss only the electromagnetic pair production by heavy ion collisions, 
the general theory can also be applied to other pair production processes in external fields, 
for example, pair creation due to bremsstrahlung |T3 . 



II. REDUCTION USING FEYNMAN BOUNDARY CONDITIONS 

The amplitude for the A^-pair creation is given with the help of the S operator above as 

Sf^ = {f\S\z). (2) 

As we want to calculate pair creation, the initial state is the vacuum state |0) and the final 
state an A^-pair state given by 

where fc, denotes the electron quantum numbers, for example, momentum and spin, and li 
the same for the positrons. The S matrix is therefore 

Sn ■■= S{ki,li ■ ■ ■ , kN,lN) = (0| dij^bkj^ ■ ■ ■ di^hk^S |0) . (4) 

If we neglect the interaction of the electrons and positrons among each other, we are 
essentially dealing with a system of independent particles. Therefore the field operator at 
any given time is linearly connected to its value at the boundary. As we can decompose the 
field operator into creation and annihilation operators, the same is also true for them. 

Using retarded boundary conditions this means that the creation and annihilation oper- 
ators for electrons and positrons at t — +oo are connected through a unitary matrix with 
those at t — >■ — oo. In the interaction picture this transformation is expressed as 

= ^Pi + ^pj ^^t^ (^^) 

j>0 j<0 

d,s = Y.<^l,sht + Y.ai,sd,. (5b) 

i>0 j<0 

The main idea is now not to use these boundary conditions but Feynman boundary 
conditions instead. These are of a mixed retarded and advanced type. With these the 
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electron creation operators for t — > +00 are connected with the electron creation opera- 
tors for t —00 and the positron annihilation operators for t +00. This corresponds 
to the Stiickelberg-Feynman interpretation of electrons moving forward, positrons moving 
backwards in time. Again in the interaction picture this is written as 0,|l^,|l^ 



q>0 q'<0 



(6) 



It has been shown already by Feynman that these boundary conditions are completely 
equivalent to the retarded ones, but that they are better suited for QED p[. For this 
relation to be applicable, we assume also that the external field vanishes asymptotically. 

Similar relations exist for all other creation and annihilation operators also; Eq. (|^) is 
the only relation we will need in the following. Applying it to we get for S^: 



Sn — (0| di^hkj^ ■ ■ ■ hk^di^ 



X 



gi>o 



g;<o 



|0). 



(7) 



The hq^ operator annihilates the vacuum state, therefore we drop this term: 

= (0| di^hk^ ■ --bkAi J2 4^ddMS\0) 



q'i<0 



= 4q[ (0| dinh^ ■ ■ ■ hkA^d^,S |0) . 
gi<o 

Using the fact that h^^ anticommutes with di^ and (i^, we get 

Sn=Y. 4,~q[ (0| dij^bk^ ■ ■ ■ di^d^ihk^S |0) . 



(8) 



(9) 



Now one replaces 6^2 using again Eq. (D and then commutes the next hk^ to the most 
right. Doing this for all hk^ we finally get: 



S 



N 



E 



q[-q'M<^ 



(10) 



The vacuum expectation value can be calculated easily using the anticommutation rela- 
tions 



di- , d\ 



6. 



h,q', 



(11) 



in order to move d"^, to the left, until it reaches the left side, where it annihilates the vacuum. 
Doing this for all creation operators we see that we get a nonvanishing result only, if every 
q[ is equal to some Ij, and the overall sign of the vacuum expectation value is equal to the 
signum of the permutation to arrange h - ■ - In in the same sequence as the corresponding 
q[--- q'j^. We get 



4 



{0\di,...di,d^,...d^,S\0)^ 



In 



(0| S |0) sgn((7)5(Z,(i), gi) • • • 5(Z,(Ar), g^v) 



(12) 



and therefore for Sj^: 

SN{ki, /i, ■ ■ ■ , k^, /jv) = 

{0\S\0)j:sgn{a)sll^^^...st;,^^^^, 

where the sum consists of A^! terms. This can also be written formally as 

(0|5|0) det 



S 



N 



(13) 



(14) 



This is the main result, which we will use. The A'^-pair creation amplitude can be written 
as an antisymmetrised product of the s^f times the vacuum amplitude. It is easy to see 
that the s^;~ are just the one-pair creation amplitudes divided by the vacuum amplitude. 
In the following we call this amplitude the "reduced" amplitude. 

Normally the vacuum amplitude does not show up in calculations, as it is often of mag- 
nitude one, for example, for time independent systems. In our case the external fields are 
explicitly time dependent, therefore the vacuum amplitude can not be neglected. This can 
be understood easily, if one notices that the vacuum amplitude is the amplitude for the 
"no pair creation process", that is, the probability amplitude for the vacuum to remain 
the vacuum. As pair creation does occur, this probability must be smaller than one and 
therefore the amplitude has to be of magnitude less than one. Also we see here that the 
vacuum amplitude appears as a factor in all A^-pair amplitudes and therefore reduces all 
these amplitudes by the same factor. 

From this we can get easily an approximate Poisson distribution for the total A^-pair 
probability, which has also been found by earlier calculations. For this we need the absolute 
value squared of S^. As Sn consists of A^"! terms, we get a total of (Nl)'^ terms. In most 
of these terms we have the product coming from two different permutations. Normally 
the produced electrons and positrons are correlated to some extent with each other. If 
we therefore assume that terms, where two different permutations appear, are much smaller 
than those with the same permutation and can therefore be neglected, we get only A^! terms: 



(15) 
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rfci'<T(i) 




s+- 



Summing or integrating now over all states we get for the total probability 

1 



total 



(N) 



m 

X E E 

ki,li,---,kN,lN 



msm 



'fcl'cr(l) 



(16) 



where the factor 1/ (A^!) has been introduced in order to correct for the multiple summation 
over electrons and positrons. As every term in the sum over a gives the same result, we get 
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Ptotai(iv)^-i(oi 5 10)1^5: 



_ 1 
~ m 

= P(0) 



ki,h 



E 

k,l 



■='kl 



k\r,lf 



N 



tdR 
-'total 



N 



(17) 



the total reduced one-pair 



with P(0) the "no pair" probabiUty and Ptotai = Y^k,i -^h 
probabihty. This is just a Poisson distribution, and P{0) can be calculated using the fact 
that the sum over all N has to be one: 



\pR 



= P(0) exp 



P. 



R 

total 



1. 



(18) 



Prom this we finally get 



-Ptotai(A^) = exp 



-p. 



R 

total 



pR 
-"^ total 



N 



(19) 



The only approximation that was necessary in order to get this result was the neglection 
of all "exchange terms" in the calculation of the probability. We show later that at least for 
the two-pair creation this neglection seems to be justified. 

Therefore for the calculation of the A^-pair creation probability it suffices to calculate the 
reduced one-pair creation probability. All other probabilities are given then by the Poisson 
distribution. If one wants to go beyond the Poisson distribution, it suffices to calculate the 
reduced one-pair creation amplitude and the vacuum amplitude. 

This derivation is rather simple, but it is not useful for explicit calculations. Especially 
the nature of the processes contributing to the reduced pair creation amplitude as well as 
the vacuum amplitude remains unexplained. Therefore we derive the same result in the next 
section using the perturbation-theory expansion of the S operator together with the Wick 
theorem. 



III. REDUCTION USING PERTURBATION THEORY 

In this section we use the field operator ^'(x) directly instead of the creation and anni- 
hilation operators as in the previous section, that is, we are going to calculate 

SN{yi,---,yN,y[,---,y'N) = {f\S\t), (20) 

where the initial state is again the vacuum |0), and the final state is 

I/) = ^^-'\yN)^^-\y'N) ■ ■ ■ ^^^\yi)^^-\y[) |o) . (21) 

Here we have made use of the frequency parts of the field operators for electrons and 
positrons. The electron coordinates are yi, the positron coordinates y^ and the time co- 
ordinate is the same for all coordinates and assumed to go to oo. S'jv gives the probability 
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amplitude in occupation number space to get electrons finally at yi and positrons at y[. We 
rewrite the final state a little bit by introducing normal ordering of the field operators: 

(/I = (0| : ^{y[)^{yi) ■ ■ ■ ^{y'^)<i>{yN) : , (22) 

where we can use the whole field operators instead of the frequency parts, as the normal 
ordering together with the vacuum state guarantees that only the appropriate part of each 
field operator contributes. Throughout the calculation in this section the Dirac indices will 
be suppressed. In the final form we arrange the operators in the appropriate form, which is 
not always possible at intermediate steps, but generally Dirac indices can be looked at, as 
if they are incorporated into the coordinates. 

Using the form of the 5* matrix from Eq. (|l]) we get 

S{yir--,yN,y'ir--,y'N) = 

{Q\:^{y[)^{yr)---^{y'^)^{yN): 

xTexpi^Jd^x : ^(a;)[-ze4(2;)]*(a;) :| |0). (23) 

As the time, when we detect electrons and positrons, is oo, the time ordering can be extended 
to the whole expression. 

We use now the Wick theorem in the following form [0, where Oj is an arbitrary 
fermionic operator: 

(0| TOi ■■ -Om |0) = E sgn(7r(j,0)(0|TO,A|0) 

X {0\TO,---Oya---OM\0). (24) 

Here i is a fixed index, '^{jii) is the permutation to put Oj and Oj to the left of all other 
operators, and the hat in the vacuum expectation value means that these operators have 
been removed from the expression. In our case the Oj are either or \l/ and we can use the 
fact that 

(0| rvl/(x)^(7/) |0) = (0| r^{x)^{y) |0) = (25) 

and 

{Q\r<^!{x)^{y)\Q) = iSF{x-y), (26) 

where Sp is the usual Feynman propagator. As we have also normal-ordered products, the 
sum over j in this case also excludes those operators, which are in the same normal-ordered 
product as the Oj. 

In a first step we apply Eq. (p4D to '^(jjn) in the expansion of Sn'- 

Im := j dxi--- dxM (0| T : ^>{y[)m{yi) ■ ■ ■ ^{y'j^)^{yN) : 

X : ^(xi) He 4(xi)] ^(xi) : ■ ■ ■ : ^(xm) [-ie 4(xm)] ^(^m) : |0) 

„ M 

= I dxi--- dxM E (o| T^iyN)^{x^) |0) He 

1=1 

X (0| r : ^{y[)^{y,) ■ ■ ■ Wr,) : ^{xi) : (1) : ■ ■ ■ : ^ : ■ ■ ■ : (M) : |0) 
= M J dzi dxi ■ ■ ■ dxM-i SpiyN — Zi)e 4(^i) 

X (0| r : ^(2/i)^(z/i) ■ ■ ■ ^Wn) ■ ^(^i) : (1) : ■ ■ ■ : (M - 1) : |0) , (27) 
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where we have used the fact that each term in the sum is identical after an even permutation 
and a renumbering of the variables Xj. This integral is just one term in the expansion of the 
exponential in S and applying this reduction to every term we get 

Sn = J dz, Sp{yr, - z,)e 4{z^) (0| T : ^{y[)^{yi) ■ ■ ■ ^{y'^) : ^{z,) 

xexpl^Jdx : ^(x) He4(x)]^(x) :| |0). (28) 

In a next step we apply the Wick theorem to "^{zi) and get similar as in the previous 
case for the integral 

:= 1 rfxi ■ ■ ■ dxM (0| T : ^{y[)^{yi) ■ ■ ■ ^{y'^) : ^(z,) : (1) : ■ ■ ■ : (M) : |0) 

N 

= -iSpizi - y[) J dxi--- dxM 

i=l 

X (0| T : ^iy[)^iyi) ■ ■ ■ ^y^^iy^) ■ ■ ■ ^(y'r,) : : (1) : ■ ■ ■ : (M) : |0) 

+M J dz2 dxi ■ ■ ■ dxM~i Sf{zi - 2:2)6 4{z2) 

X (0| T : ^{y[)^{yi) ■ ■ ■ ^{y'^) : ^(z^) : (1) : ■ ■ ■ : (M - 1) : |0) , (29) 
which is again a term in the expansion of S^, so that our amplitude is 

N 



Sn = ^-i j dzi SpiyN - zi)e 4{zi)SFizi - y[) 

i=l 

X (0| T : ^{y[)^{yi) ■ ■ ■ ^ly[)^{yi) ■ ■ ■ ^{y'j,) : 
xexpl^Jdx : ^'(x) [-ie4(a;)] *(a;) :| |0) 

+ J dzidz2 SpiyN - zi)e 4izi)SFizi - Z2)e 4(2:2) 

X (0| T : ^iy[)^iyi) ■ ■ -^{y'^) : ^(^2) exp [J dx : ^(x) \-ie 4(x)] ^(x) :| |0) . (30) 

The vacuum expectation value in the first term is just the (A^— l)-pair creation amplitude, 
that is, 

= / '^^1 SpiyN - zi)e 4-{.zi)Sf{zi - y[) 

i=i 

xsgn(a)S7v-i(yi, • ■ ■ , yN-i, y[r--,yir--, v'n) 
+ J dzidz2 SpiyN - zi)e 4(^1)5*^(^1 - Z2)e 4(^2) 

X (0| T : ^{y[)^{yi) ■ ■ ■ ^{y'r,) : ^(^2) exp S^J dx : f (x) [-ze 4(x)] ^(x) :| |0) , (31) 
where a is the permutation to rearrange the ^{y'j) and "^{yj) to the standard form. We 



can use Eq. (pOf ) again for the vacuum expectation value in the second term. Doing this 



recursively and assuming that the series converges we finally get 
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N 

Sn = Y^sgn{a)SN-i{yi, - ■ ■ ,yN-i,y[, - ■ ■ ,y'i, - ■ ■ ,yN) 

i=l 

X (-i) y dzi SpiyN - zi)e 4{zi)Sf{zi - y[) 

+ j dzidz2 SpiyN - zi)e 4{zi)Sf{zi - Z2)e 4{z2)Sf{z2 - y[) 

+ •••]■ (32) 

The infinite series is just the normal perturbation theory expansion for the pair creation 
without any disconnected parts. We identify it with the reduced one-pair creation amphtude 
S^{yN, y'j). Therefore we have found a way to reduce the iV-pair amplitude to a (A^ — l)-pair 
amplitude. Recursively we can write for Sj^: 

SN{yi,---,yN,y[,---,y'N) = 

SoY.mn{a)S''{y^,y',^,^) ■ ■ ■ S^{yr,,y',^r,)) (33) 
with the vacuum amplitude 

5*0 = (0|rexp|y dx : ^ (x)[-ie 4{x)]^ {x) :| |0) 

= {0\S\0). (34) 

This is identical to the result of the last section. Therefore we have found an expression 
for S^, that is, for . It is just the perturbation theory expansion of a fermion interacting 
an arbitrary number of times with an external field without any disconnected parts. 

Finally we want to find also a perturbation theory expression for 5*0. Expanding the 
exponential we get 

OO -1 

^0 = E TH^M, (35) 

M=0 -'^■^ • 

with Im defined as 

•= J dxi ■ ■ ■ dxM 

X (0| r : ^{xi) [-ie 4{xi)] ^(xi) : • • • : (M) : |0) . (36) 

We are now using the Wick theorem again in order to reduce this expression. Applying 
it to ^{xi) we get 

Im — {—i) {M ~ ^) j dxidx2 ■ ■ ■ dxM 
xe 4(xi)5'f(xi - X2)e 4{^2) 

X (0| r^(xi)*(x2) : (3) : • • • : (M) : |0) . (37) 
Similar we get for the following expression 
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Jm{xi, X2) := jdxs--- dxM (0| T^{xi)^{x2) : ^(^3) [-ie ^(xg)] ^{x^) : • • • : (M) : |0) 
= -tSF{x2 - xi) Jdxs--- dxM (0| r : (3) : • • • : (M) : |0) 
+(M ~ 2) y dxsdx4 • ■ ■ dxM iSF{x2 — x^) [—ie ^(^3)] 

X (0| T^{xr)^{x^) : (4) : ■ ■ ■ : (M) : |0) . (38) 

That is, 

iM^-i{M-l) j dxidx2 e 4-{xi)Sf{xi- X2)e 4-{x2)Jm{xi,X2) (39) 



and 



Jm{xi, X2) = -iSF{x2 - xi)Im-2 + (M - 2) j dxs Sf{x2 - x^)e 4-{x^)Jm-i{xi, x^). (40) 

Using both relations in order to reduce Im in terms of 1^ we get 

iM = {M- 1) C2/m_2 + (M - 1)(M - 2) C3/M-3 
+ • • • + (M - 1) • ■ • 2 Cm-Ji 

+(M-l)---2 1 CmIo, (41) 



with Cjv defined as 



Cm : 



N '■= I, — 1) y <ia;i • • • dxN 



Tr 



e 4(xi)5'i^(a;i - ^2)6 4(^2) 



xe ^{xn)Sf{xn — Xi) 



(42) 



writing explicitly the Dirac indices to get the trace. For ease of writing we define also 
Co — Ci — and use that h — and Iq — 1- With these we can write more compactly 

^Ji^^^'W~ky.- ^^^^ 

This reduction formula can now be used in order to find an expression for So- For this we 
write Im in a series ordered by the number of Q's (that is, as we will see later, the number 
of loops) that appear in each term after a complete reduction. We see that the formula 
above just gives us a recipe, how the next order can be derived from the lower one. Only 
one term with no in it exists: 

= 1, (44) 

giving therefore in zeroth order for 5*0: 
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0(0) 
•^0 



(45) 



Using this only nonzero in Eq. (^3D we get in first order only a contribution if /c = M 



?(i) 
'a/ 



M! 



and for Sn 



In second order we get 



"^0 



(1) 



E 

M=0 



M! 



E 



Cm 



(46) 



(47) 



?(2) 



M! 



Af 



M ^ 



?(i) 

'Af-A: 



fe=l 



{M-k)\ 



~M ~k M -k 



M 2 



M! 1 ^ ^ 

fc=i 

E 



+{M-k) ^^-^ 



k M-k 



M-k k 



Ml ^ Ck CM-k 



k=l 



k M-k' 



(48) 



where we have changed the summation index from k to M — k in one step. For 5*0 we get 



C(2) 



1 ^^Ck Ca 



oo 1 

Af=0 ' ^ Af=0 fc=l 



oo 



E 

.fc=i 



oo 



E 

./=1 



5, 



(1) 



(49) 



Using the same transformations we can express also all higher terms of 5*0, for which we get 



nl 



5, 



(1) 



(50) 



so that we finally get the vacuum amplitude as 



5, 



(1) 



exp 



oo 

E TT / dxi ■ ■ ■ dXM 



M=2 



M 



So = exp 

xTr [e 4{xi)Sf{xi - X2) ■ ■ ■ e ^{xm)Sf{xm - Xi)]^. 



(51) 



This is a well known result [^-[Tl|,[l^. Its interpretation is straightforward. The sum in the 
exponential is just the sum over all single loops, whereas the exponential accounts for the 
fact that we can have also two or more loops. 

With this we have found now perturbation-theory expressions for all terms occurring in 
the reduction formula Eq. (P^D, that is, for the reduced one-pair amplitude and the vacuum 
amplitude. 

A final remark has to be made about the vacuum amplitude Sq. Already Schwinger 
found out that the expression for Sq is not finite due to infinities in the imaginary part of 
S'q^ hi (see also ||12|). But he also mentioned that in the total probabilities we need only the 
absolute value squared of 5*0, where only the real part of S^^ contributes, which is finite. 
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IV. COMPARISON WITH EARLIER RESULTS 



The results of the previous two sections suggest the following picture of the pair pro- 
duction: In the Feynman picture the A^-pair production can be described by two forms of 
processes (Fig. |lD. positron lines enter the interaction region coming from the future. 
They interact with the external field an arbitrary number of times, where they may change 
also their direction in time. Finally they leave the interaction region as electron lines mov- 
ing into the normal time direction. Besides these processes, which are characterized each 
by a continuous line coming from and leaving to the future, there are also processes, which 
consists of closed loops. As they remain entirely inside the interaction region, they are not 
visible as physical processes. These closed loops form the vacuum amplitude. From this 
picture it is also clear that the vacuum amplitude is the same for all A^-pair amplitudes 
and therefore a common factor in all of them. The fact that we are dealing with multiple 
particles only shows up in the calculation through the antisymmetrisation with respect to 
all electrons (or equivalent to all positrons). 

This means that for a complete calculation of the A^-pair creation probabilities we need 
to know only the reduced one-pair creation amplitude together with the vacuum amplitude. 
Both can be calculated in principle using perturbation theory, where the reduced one-pair 
amplitude is identical to the usual perturbation series result and setting the vacuum ampli- 
tude to one, that is, neglecting all diagrams with disconnected parts. 

If one neglects the exchange terms in the calculation of the probability, it suffices to know 
the reduced one-pair creation probability, as this is the only result needed for the Poisson 
distribution of the multiple-pair probabilities. 

Note that according to the rules of the Feynman diagrams we have neglected the antisym- 
metrisation and therefore the Pauli principle for all intermediate states. There is a deeper 
reason for doing so, which has to do with the connection of the higher-order processes and 
the occurrence of the vacuum processes through antisymmetrisation. This is, for example, 
discussed in detail in [0,0 . Let us look, for example, at a typical higher-order diagram in 



perturbation theory, where we would expect corrections because of the Pauli principle. In 
Fig. we would expect a deviation as we have two electrons, which are not allowed to be in 
the same state. (The same is also true for the two positrons.) But it can be shown that the 
contributions from this process are just canceled by those of the vacuum correction shown 
in Fig. 1^, where the electron in the loop and the produced electron are not allowed to be in 
the same state also. Both processes are connected with each other, as we get one from the 
other, if we exchange the two electron lines or the two positron lines. This result is of a very 
general nature, so that the antisymmetrisation of all intermediate particles can be dropped, 
if we also include the vacuum processes (again without antisymmetrisation). 

This shows that the higher-order multiple-particle corrections to the one-pair creation 
are connected with the vacuum process. Therefore we think that these processes should be 
studied in more detail. Especially as we know that the vacuum processes are not negligible, 
as they decrease the reduced probabilities larger than one to values less than one. 

The inclusion of these higher-order processes is straightforward in our model. As the 
Poisson distribution needs only the neglection of the exchange terms, we can include the 
higher-order terms into the calculation of the reduced probability and use this result then in 
the Poisson distribution for the multiple-pair probabilities. This procedure is consistent, as 
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it includes automatically the higher-order corrections to the vacuum amplitude. Therefore 
the calculation changes only the reduced probability used in the Poisson distribution but 
not the explicit form. 

Let us compare this result with the other models AH three models get a Poisson 

distribution for the multiple-pair creation as well but only based on the summation of a 
restricted class of diagrams. Studying further these approximations one finds that all of them 
are essentially "quasi boson" approximations. This means that they have as fundamental 
processes a pair creation and a pair annihilation process (neglecting at the moment the 
Coulomb scattering). Combining these processes all previous models assume that only pairs, 
which have been produced together in a creation process, can be annihilated. Therefore the 
electron and the positron are seen as an unbreakable pair, which behaves more or less like a 
boson. And for bosons in an external field the Poisson distribution comes out exactly |p!0 



Even the inclusion of Coulomb scattering in some of these models does not change the "quasi 
bosonic" nature of them. Our calculation indicates that there are higher-order processes in 
which electrons and positrons do not behave as "quasi bosons" . These processes are multiple- 
particle processes in the sense that more than one electron or positron are essentially needed 
in an intermediate step. 

The fact that these are multiple-particle effects also shows the advantage of using the 
Feynman boundary conditions instead of the retarded boundary conditions used in the Dirac 
sea picture. In the Feynman picture the production of a pair is described by a continuous 
fermion line, where the inclusion of multiple-particle effects presents no difficulties. The 
advantage of the Dirac sea picture is mainly that its particle hole interpretation allows to 
treat the creation of one electron-positron pair in a single particle formalism. As the higher- 
order processes need the existence of more than one electron and one hole, this advantage 
of the Dirac sea picture is no longer existent then. 

The other advantage of the Feynman boundary conditions is that the different processes 
can be separated from each other. Especially the vacuum processes are independent of the 
reduced pair creation processes. This is not the case in the Dirac sea picture, where vacuum 
processes and pair creation processes can not be separated in the intermediate steps but 
only in the final state. 

In the following we show that the processes of the type of Fig. |^ are the lowest-order 
corrections to the one-pair creation processes if we neglect Coulomb scattering terms. For 
this we restrict ourself to second-order Magnus theory. 

V. THE S OPERATOR IN MAGNUS THEORY 

The Magnus theory can be seen as an expansion in the interaction time. For large values 
of 7 the interaction time of the two heavy-ion fields is very short, therefore the use of this 
approximation seems to be justified. Up to second order the S operator is given by [0,0 



S = T exp 
^ exp 



/ + 00 
Hj{t) dt 
-oo 

/+00 1 r+oo rt2 

Hi{t)dt + -i-tf / dh / dt, [Hj{t2), Hj{ti)] + ■■■ . (52) 
-OO ^ J — OO J — CO 
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The commutator of the if/'s for ti < t2 is the difference of the time-ordered and the 
anti-time-ordered product. We rewrite S as 

S = expl -i / Hj{t)dt + -{-if / dt2 / dti 

I J —OO 4 J — QO J — OO 

xT[Hi{h)Hi{t^)] - A[Hj{t2)Hj{h)] + ■ • j (53) 

= expjy d'^x e : ^{x)[-ie 4(x)]^(a;) : 

d^xi d\2 {T-A): ^>{x2)[-ie^{x2)]^{x2) : : ^ {xi)[-ie ^{xi)]^ {xi) :|. (54) 

Now we use again the Wick theorem in order to put the field operators into normal- 
ordered form. We use it for the time-ordered and anti-time-ordered products in the form 

TAB = :AB:+ {Q\TAB\Q) = : AB : + {TAB) , (55a) 
AAB = :AB:+ {<d\AAB\Q) = : AB : + {AAB) . (55b) 

As the first term in S is in normal-ordered form already, only the second one has to be 
rearranged. For the time-ordered product we get 



+ : ^2^1 ■■ {t^2^i) + (Tf 2^1) {T^2^i) , (56) 
and the same for the anti-time-ordered product by replacing T with A. We get for S 

S = expjy" [-ie 4{x)a^{x)] : ^a'^p ■ +^ J dxidx2 [-ie 4{x2)ap][-ie 4{xi)./s] 
x[: ^p{x2)^^{xi) : (V^a{x2)^5{xi))+ : ^a{x2)^5{xi) : {V^p{x2)^-y{xi)) 

+ (T^^{x2)'^5{Xi)) (T^p{x2)^-y{Xi)) - (A^^{x2)'^5{Xi)) (^^^(x2)^'^(Xi))] |, 

(57) 

where we have defined P as T — The last two terms can be rewritten to give 

[V^^ix2)^5ixi)) (T^pix2)^^ixi)) + (A^^ix2)^sixi)) (v^ f,ix2)^ ^(xi)) . (58) 



Exchanging the integration of dxi and dx2 we see that the second and the third term of 
Eq. (pTl) can be combined and also the last two terms, so that we get 



5" = expjy [-ie 4{x)ai3{x)] : \l/a^/3 : +^ J dxidx2 [-ie 4{x2)af3][-ie 4{xi)-f5] 
X 2 : <ilp{x2)^jixi) : (V^a{x2)'if5ixi] 



+ ((T + A)^pix2)^^{xi)) {V^^{x2)^5ixi))^ j. (59) 

Now we transform S into momentum space. For this we need the Fourier transform of 
the field operator as well as of the vacuum expectation values. 
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The field operator and its conjugate can be written as 



b{p, s)u{p, s) exp{—ipx) 



+d^{p, s)v{p, s) exp{ipx) 



b'^{p, s)u{p, s) exp{ipx) 



+d{p, s)v{p, s) exp{—ipx) 

where we have introduced dp, which is defined as 

d^p 



dp :- 



(27r)3/2(2po)^/^ 



(60a) 



(60b) 



(61) 



In the following we need also the factor of dp alone together with the Lorentz invariant phase 
space, therefore we define also 



7(P) 



1 



and 



dT{p) : 



(27r)3/2(2po)i/2^ 
d^p 



(27r)32po' 



(62) 



(63) 



With the usual relations between vacuum expectation values, propagators, and singular 
function, one gets for the vacuum expectation values ]TP|-pi|] 



(0| T^(x)^(a;') |0) = iSpix - x' 
d*p 1^ + m 



(2-7?)^ p2 _ j7^2 _|_ 



exp(— ^^(x — x')). 



(64a) 



(0| ^^(x)^(x') |0) = -iSa{,x - x' 
d^p ^ + m 



—I 



Combining both we get for T) and T + A: 



exp(— ^^(x — x')). 



(0| V^{x)-^{x') |0) = i [Sf{x - x') + Sa{x - x') 



J (27r)* \p^ — m 



+ It p^ — — le 



X exp(— 2p(x — x')) 



2i / ,^^.. {:^ + m) J''^' exp{-ip{x-x')), 



p2 _ ^2 



(64b) 



(65) 



where P.P. denotes the principal part of the integral, and 
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(0| (T + ^)^(a;)^(x') |0) = i [Spix - x') - SAix - x')] 
dT{p) + m) exp{—ip{x — x')) 

dT(p) — m) exp(zp(x — x')), (66) 

which describes the propagation of on-shell electrons and positrons. 

These forms of the vacuum expectation values are put into S together with the decom- 
position of the field operators and ^{x) and using also the Fourier transform of the 
external field 

A{x) = I A(g) exp(— zgx). (67) 



; {2'nY 

The integration over the coordinate space can be done, giving deltafunctions for the 
momenta. Reducing the momentum integrals with the help of these, we finally get 



S = expj— J dpi dp2 



X : h^{pi)h{p2) : u{pi) 4{pi - P2)u{p2) 

+ : b^{pi)d^{p2) : u{pi) 4{pi + P2)v{p2) 
+ : d{pi)h{p2) : v{pi) 4{-Pi - P2)u{p2) 

+ : d{pi)d'^{p2) : v{pi) 4{-Pi +P2)v{p2) 
dp 



—ie^ J dpi dp2 



P.P. 



: 6+(pi)6(p2) : u{pi) 4{Pi - p) + 

p.p. 

+ : b^{pi)d^{p2) : u{pi) 4{pi - p) (^ + m) — 

pZ _ ^z 



4{P - P2)U{P2] 



+ : d{pi)h{p2) ■■ v{pi) /(-pi -p){^+ m) 



P.P. 



p2 _ ^2 



+ : d{pi)d'^{p2) ■■ v{pi) 4i-Pi - + fn) 



+i- 



-'J 



dp 
(2^ 

dp 
(2^ 



dr{p') Tr 
dT{p') Tr 



m) j^{p — 



P.P. 

p2 _ ^2 



+ m 



4{P+P2)V{P2) 

4{P-P2)u{p2) 

4{P + P2)V{P2) 
P P 



p2 _ ^2 



P.P. 



(/ - m) -p)(^ + m) 

pz _ ^z 



4{p + p') 



(68) 



Again, we have not written explicitly the spins of the leptons, which we look at as 
included in the momenta. 

The interpretation of the individual terms is straightforward (see Fig. Terms of the 
form : b'^b : are electron Coulomb scattering terms, : dd~^ : the corresponding positron 
scattering terms. : b'^d'^ : corresponds to pair creation and : db : to pair annihilation. The 
last two terms have no operators in them. They correspond to the lowest-order vacuum 
corrections. As they are only numbers, they commute with all other terms. 
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VI. APPLICATION TO ELECTROMAGNETIC PAIR CREATION 



We apply this result, which is true for an arbitrary external field, to the special case of 
pair production in heavy-ion collisions. In this case there are some kinematical restrictions, 
so that some of the terms in S can be dropped. Especially there is no pair creation or pair 
annihilation in first-order Magnus theory and we can replace the principal value integrals 
in the pair creation and annihilation term in the second order by ordinary integrals, as the 
intermediate state is not allowed to be on-shcU in this case. 

As we want to concentrate on the multiple-particle effects we will neglect also all Coulomb 
scattering terms. It is generally assumed that these terms do change the differential prob- 
abilities, but are only of minor importance for the total probabilities. Also their infiuence 
seems to become smaller for higher energies. Therefore we will neglect all terms of the form 
: b~^b : and : dd'^ :. Also the vacuum terms will be dropped, as we are only interested in 
the calculation of the reduced amplitude, where vacuum corrections do not appear. The S 
operator with these approximations is 



The integrals over p in this expression are just the pair creation and pair annihilation 
amplitude in second-order Born approximation. We define the pair creation and annihilation 
potential V and V* as 




X 6+(pi)d+(p2) 

Xu{pi) 4{pi-p)— 4{P + P2)V{P2) 

+d{pi)h{p2) 




(69) 




Xu{pi)4.{pi -p)^ o 4{P + P2)V{P2) 



(70a) 




xv{pi)4{-pi -p)^— — ^4(P-P2)^i(p2) 



(70b) 



(where V* is just the complex conjugate of V) in order to write S as: 




(71) 
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The matrix element of the pair creation in second order is connected with the potential V 
through 



Mip,q)^-iVip,q). (72) 

For the reduced one-pair creation amplitude we have to calculate 

{f\S\i)^{0\diqf)bipf)S\0). (73) 

We are now expanding the exponential in S in order to get the different contributions to 
the one-pair creation. In lowest order we expect to get back the second-order Born result. 
We get 

(/I \i) ^-ijdp dq V{p, q) (0| d{qf)b{pf)b+{p)d+{q) |0) 
= -h{Pfh{Qf)V{Pf, Qf) 

^M(pf,qf)-f(pf)^{qf), (74) 
The differential probability in first order therefore is 

P{pf,qf) = \M{pf,qf)f^\pf)^\qf) (75) 
and the total probability is given by 

Ptotai = / \M{pf, qf)\' dr{pf)dr{qf), (76) 

which is indeed the Born result. 

As the vacuum expectation values, we get from higher orders, vanish, if the number of 
creation and annihilation operators is not equal for each kind of particles, we see easily that 
only odd orders in the expansion of S contribute. The next order is therefore the third one, 
where we get 

(/I S^^'^ \i) = ^ j dpi dp2 dp3 dqi dq2 dq^ 

X V{pi, qi)V*{p2, q2)V{ps, q^) (0| d{qf)b{pf)b+{pi)d+{qMQ2)b{p2)b+{ps)d+{qs) |0) 

+V*{pi, qi)V{p2, q2)V{ps, qs) (0| d{qf)b{pf)d{qi)b{pi)b+{p2)d+{q2)b+{ps)d+{qs) |0) 

(77) 

Here we have used the fact that only two of the eight possible combinations do not vanish 
trivially. Again the vacuum expectation value can be calculated using the anticommutation 
relations of the particle operators. The first expectation value is 

{0\d{qf)b{pf)b+{pi)d+{qi)d{q2)b{p2)b+{ps)d+{qs) |0) = 

^{Pf -Pi)S{P2-P3)S{qf - qi)S{q2 - qs) (78) 

and the second one 
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(79) 



Terms of the form 6{p2 — P3)5(q'2 — Q's) describe closed fermion loops, and we drop them in 
the calculation of the reduced amplitude. Therefore only two terms remain. In third order 

we get 

(/I S^'^ \t) = -47(P/)7(g/) / dTip^)dT{q^)Vipf,q,)V*ipuqi)Vipuqf). (80) 



3! 

Combining both results we get up to third order 

(/I 5(1+3) \i)=^ipf)^iqj) ^M{pf,qf) + ^ j dV{p,)dV{q,)M{pf,q,)M*{p,,q,)M{p,,qj) 

(81) 

The interpretation of this higher-order process is straightforward. If we keep in mind 
that M{p,q) is the amplitude for the production of a pair and M*{p,q) the amplitude for 
the annihilation of a pair, we get the Feynman diagram as shown in Fig. |^. These are just 
those higher-order processes, we already discussed before in Sec. |I^ (see Fig. |^). Two pairs 
are created and then one of the electrons annihilates with the other positron, so that finally 
only one electron and one positron remain. The Magnus theory restricts these intermediate 
particles to the mass shell, so that we only need the on-shell M, whereas they may be 
off-shell in the general case. 

Using the same technique higher orders can also be calculated. The explicit calculation 
of the fifth and seventh order gives 

16 r 
{f\S^''^\^) =-j{pfUqf) / dT{p,)dT{q^)dT{p2)dT{q2) 



5! 

xM{pf, qi)M*{p^, gi)M(pi, g2)M*(p2, g2)M(p2, g/) (82a) 
272 r 

{f\S^'^\€) = —j(pj)j(qj) J dT{p,)dT{q,)dT{p2)dT{q2)dT{p,)dT{qs) 

xM{pf, q,)M*ip,, qi)M{p,, q2)M*{p2, q2)M{p2, qz)M*{p^, q3)M{ps, qf), (82b) 

where in the calculation of the vacuum expectation values also more complicated fermion 
loops occur, which have to be neglected. The following properties are remarkable, as they 
are true for all higher orders: First all contributions coming from the same order are of the 
same type and can therefore be summed. Second all amplitudes have the same sign, that is, 
they are all added coherently. Therefore they increase all the reduced total probability, no 
cancellation or reduction occurs. This same sign can also be understood from the connection 
of the higher-order processes with the vacuum corrections. Vacuum loops contribute with a 
negative sign to the amplitude (see Sec. fT| and JI^)- As our diagrams are connected with 



a vacuum diagram by the permutation of some lines, it is clear that it has to have the same 
sign as the lowest-order diagram. The general form of the nth order therefore is 
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= c„7(p/)7(g/) J rfr(pi)dr(gi) ■ ■ ■dr(p(„_i)/2)rfr(g(„_i)/2) 

xM{pf, gi)M*(pi, gi)M(pi, gs) ■ ■ ■ M*(p(„_i)/2), g(„_i)/2)M(p(„_i)/2, g/), (83) 

For higher orders the number of terms occurring in the reduction of the vacuum expec- 
tation values gets rather large and the number of diagrams increases rapidly. However we 
show in the appendix how these total number of processes and therefore c„ can be calculated 
using combinatorial arguments and a recursive formula. 

Finally we calculate the lowest-order correction to the reduced total probability. Squaring 
the pair production amplitude up to third order (Eq. (|81|)) we get 



= M*ipf,qj)- 



dTip)dT{q)M*{pf, q)M{p, q)M*{p, qj] 



M{pf,qf) + - I dr{p')dr{q')M{pf,q')M*{p',q')M{p',qf] 



|Mb/,g/)r + 



+M*{pf, qf)M{pf, q)M*{p, q)M{p, qj) 
2 



dT{p)dT(q)M{pf, qf)M*{pf, q)M{p, q)M*{p, qf) 



\M{pf,qj)\' + -JdT{p)dT{q) 
xRe [M{pf, qf)M*{pf, q)M{p, q)M*{p, qj)] + ■ ■ ]l\pfh\qf). 



(84) 

Integrating over pf and g/ we get the total probability (Fig. ^ 

Ptotai = / dTipf)dT{qf)\Mipf,qf)\' + ^ I dTipf)dTiqj)dT{p)dTiq) 

xRe [M{pf, qf)M*{pf, q)M{p, q)M*{p, qj)] (85) 
=. p(B) + p(Af)_ (86) 

This result is easily generalized. The general series is of the form 

Ptotai = JdT{pi)dT{qi)\M{puqi)\^ 

+4 / dT{pi)dT{q,)dT{p2)dTiq2)Re [M{p,, qi)M*{p,, g2)M(p2, g2)M*(p2, Qi)] 
+ ■■■ 



+d2i-i / dT{pi) ■ ■ ■ dT{qi)Re [M(pi, gi)M*(pi, q^) ■ ■ ■ M{pi, qi)M*{pi, q,)] 



+ 



(87) 



The coefficients dn are again calculated in the appendix. They can be derived easily from 
the coefficients c„ of the amplitudes. 
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VII. CALCULATION FOR IMPACT PARAMETER ZERO 



For the calculation of the lowest-order multiple-particle correction to the one-pair cre- 
ation probability we make use of the analytic form of the matrixelement for impact parameter 
b zero. Details of this calculation can be found in a previous publication 0, therefore we will 
only review rather briefly the main properties. The Feynman diagrams contributing to the 
lowest order Born calculation are shown in Fig. ^ where (1) and (2) denotes the interaction 
with the external field of ion 1 and 2 , respectively. The matrix element is 

M{p, q) = —ie^u{p) 



X 



+ 



(2vr)4 



J^2 _ ^2 



The external fields are given by 



^2 _ j^2 



f(g). 



4^'^)(g) = -2vrZ(i,,)e<'^)5(gn(^'^)) 



g2 



(88) 



(89) 



u^^'"^^ are the four velocities of the ions and F{q^) is their form factor. All calculations are 
done for collisions of the same type of ions and in the center of mass system. For the form 
factor a simple dipole form factor of the form 



dipole 



A2 



q^ 



(90) 



with A = 83MeV, has been used together with a second form factor, which is a linear 
combination of two dipole form factors. In the previous publication the integral has been 
solved in order to calculate the lowest-order Born result. Here we use this analytic result 
for M in Eq. (p6|). The usual way to calculate expressions of this kind is to rewrite the 
spin summation as a trace over gamma matrices. But as each of the M consists of two 
diagrams, we finally get a total of sixteen different traces, which all give large expressions, 
which are not manageable by an algebraic calculation program and therefore can not be 
used for numerical calculations. Therefore we prefer to calculate M directly and to do the 
spin summation numerically. 

Several methods have been proposed in the literature how to calculate amplitudes directly 
instead of their squares p^-p^. We use the method described by Fearing and Silbar ||25|| . 
It consists of multiplying and dividing M with v{q)u{p) in order to get 



M = u{p)Mv{q) 
1 



Tr< v(q)v{q)u{p)u{p)M 



v{q)u{p) 

^^Tr{(^-m)l(l + 75A,/. 



X (^ + m)-(l+75Ap /p)M 
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Tr< m) 



4S{p,q) 

X (1 + TsAg + TsAp - XpXq /q /p) 

x(^ + m)Ml, (91) 



where Sp, Sq are the spinvectors of the electron and the positron, and Ap, are the eigenval- 
ues of the spinors with respect to the spinvector, that is, they have values of ±1. S{p, q) is 
a complex number, which has been calculated using an explicit form for the spinors for the 
standard form of the gamma matrices. In order to make the numerical calculations easier 
we have used polarisation vectors Sp and Sg, which are longitudinal vectors: 

(p.,0,0,po) (92) 



and similar for Sq. The calculation of the trace has been done with the help of the algebraical 



calculation program FORM |£6[ . All sixteen different spin combinations have been calculated 
and summed. The final expression has to be real, which has been used to test the program 
(As the final expression can be expressed as a trace, which does not contain any 75 matrices, it 
must be equivalent to an expression containing only real numbers and real scalar products, 
therefore it has to be real). The integration over four particles gives a 12- dimensional 
integral, of which one angular integration is trivial. It was calculated using a Monte Carlo 
(MC) integration routine (VEGAS [^3,^]). Together with the total probability we have 
calculated also single differential probabilities by sorting the points into appropriate bins. 
The error of the total probability is always less then 1%. An explicit error estimate is given 
for the differential probabilities. As our calculation is symmetric with respect to all four 
momenta and all four electron-positron combinations, this error has been calculated from 
the standard deviation of the four results. 

Finally we want to point out that the expression in Eq. ( ^61) can be interpreted also in 
a different way: If we calculate the reduced total two-pair creation probability exactly in 
lowest order according to the previous section, we get 

All = / dT{p{)dT{p,)dT{q{)dT{q^) 

X [M*(pi,gi)M*(p2,g2) -M*(pi,g2)M*(p2,gi)] 
X [M(pi, gi)M(p2, q2) - M{pi, q2)M{p2, qi)] 

I [ dT{p,)dT{p2)dT{qi)dT{q2) 



2 

x||M(pi,gi)HM(pi,gi)| 
-Re [M*(pi, gi)M(p2, qi)M*{p2, g2)M(pi, q^)] 

— Aotal ~ Aotal- (93) 

The first term is the direct part, which is used in order to get the Poisson distribution, 
whereas the second term is the exchange part, which we neglect in the Poisson distribution 
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(see also Fig. ||). Comparing the second expression with that in Eq. ( ^6]) we see that they are 
identical. Comparing Fig. ^ with Fig. ^ we see that we only have to change the interpretation 
of the diagram by cutting it at a different point. Therefore the result of the calculation for 
the total probability can also be used to test, whether the neglection of the exchange term 
in order to get the Poisson distribution is justified. 

VIII. RESULTS AND CONCLUSIONS 

Figure ^ shows the results of our calculations for the reduced probability as a function 
of 7. In all our diagrams we set Za = 1, as this is a common factor in all results. Also 
shown is the result for the corresponding probability of the lowest-order calculation (see 
1^). In order to test the dependence on the form factor the calculations have been done 
for the dipole form factor and the double dipole form factor. Both agree within the error 
interval, therefore we show here only the results for the double dipole form factor. We see 
that the correction is similar in size as the lowest-order result for large values of 7. Using 
realistic values for Za its importance is reduced, as the higher-order correction has to be 
multiphed by {ZaY, the lowest-order result by {Za^. In Table | we give predictions for the 
contribution of the higher-order correction to the reduced total probability. For 7 and Z 
we have used typical values for relativistic heavy-ion colliders. The higher-order correction 
increases the reduced probability for Pb and U up to about 5-10%. Therefore they should 
be observable in principle. Based on this calculation we expect that higher-order corrections 
(5th order and more) are again only of about 5-10% of the third-order results and are 
therefore corrections of less than 1% to the total probability. 

Also shown are the results for the differential probabilities as a function of the energy 
E (Fig. p!0| ) and the angle with the beam axis 6 (Fig. |1T]). In these as in all other single 
differential probabilities as well the correction follows more or less the lowest-order result. 

In Fig. [1^ we compare both terms contributing to the total two-pair production prob- 
ability. Here the first and the second diagrams of Fig. ^ are multiplied by the same factor 
{ZaY, therefore the ratio of both curves gives directly the contribution of the exchange term 
to the first term, which is used in the Poisson distribution alone. We see that the exchange 
diagram contributes only with about 1% to the reduced total probability. Also its impor- 
tance gets smaller for larger values of 7, therefore the use of the Poisson distribution seems 
to be justified. This could be different, of course, if one looks at differential probabilities in 
some region of the phase space. 

Let us summarize our results: We have shown that the A^-pair creation amplitude can 
be reduced exactly to the reduced one-pair creation amplitudes and the vacuum amphtude. 
Therefore for the calculation of all A^-pair probabilities, it suffices to calculate those two 
amplitudes. The use of the external field approximation and the neglection of the interaction 
among electrons and positrons is essential in order to get this result. We have shown that we 
get a Poisson distribution for the total probabilities very generally, if we neglect all exchange 
terms. A calculation of the reduced one-pair creation probability then suffices alone, as all 
higher-order processes are given by the Poisson distribution. 

We have found perturbation theory expressions for the reduced amplitude as well as 
for the vacuum amplitude, so that both can in principal be calculated using perturbation 
theory. None of them is of a principal nonpertubative character. 
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We have compared this result with earher calculation, where we have pointed out that 
they have mainly made use of a quasi boson approximation of the electron-positron pairs in 
order to get the Poisson distribution. On the other hand our form of the Poisson distribution 
needs only the neglection of the exchange terms, therefore higher-order processes beyond the 
quasi-boson model can easily be incorporated. 

Based on the general form of the S operator in second-order Magnus theory we have 
calculated the lowest-order contribution of such higher-order multiple-particle effects to the 
total probabilities for collisions with impact parameter b zero. The contribution were found 
to be on the 1 - 10% level for high Z, and should therefore be incorporated into the calcu- 
lations. 

Finally we compared the deviation of the two-pair creation probability from the Poisson 
result. This deviation was found to be small, therefore the use of the Poisson distribution 
seems to be justified. 



As shown in Sec. all higher-order processes in the Magnus theory without Coulomb 
rescattering terms are of the same type. The number of terms in each order is given by the 
reduction of the vacuum expectation values. The explicit reduction is not useful because 
the complexity of the expressions increases quickly. Therefore we show here how these can 
be calculated more easily. 

This is a combinatorial problem, which can be formulated in the following way: In iVth 
order we have interaction points, each of which either corresponds to a pair creation or 
pair annihilation process. In terms of the fermion lines this means that a line entering one 
point from the left also leaves it to the left and the same from the right. We need the total 
number of paths through these N points using the above condition, where the line initially 
enters from the right and finally leaves to the right. As we look only at reduced amplitudes, 
no closed loops are allowed. 

Let us find a recursive formula for the number of paths: We define A{n, i) as the number 
of distinct paths for i fermion lines (coming from and leaving to the right) going through 
n interaction points. If we add now another interaction point to the right, there are two 
possibilities: If it is a pair creation process, the number of lines is increased by one. If we 
have a pair annihilation process, we have to connect it with one outgoing and one ingoing 
line, therefore reducing their number by one. For this we have to choose one of the i 
outgoing lines and one of the ingoing lines, where we have to be careful not to make a closed 
loop; therefore only i — 1 oi them are allowed. This means that there are i{i — 1) possible 
combinations. This gives us the recursion relation 



APPENDIX A: COEFFICIENTS IN MAGNUS THEORY 



A{n + l,i) = A{n,i 



+ l)A{n,i + 1) 



(Al) 



together with the boundary conditions 



A(l,l) = l 
A{l,i) = 
A{n,0) = 



for all n. 



for i > 1 



(A2a) 
(A2b) 
(A2c) 
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With this recursion relation we can calculate A{n, 1) for all n. This is just the number of 
diagrams that we were looking for. Dividing it by n\ gives the coefficient c„ in the expansion 
of the amplitude. These coefficients are given in Table ^ For the coefficients in the reduced 
total probability we have to multiply q and c„+i_; and sum over all / e . . . n. These are 
the coefficients dn of Eq. (|87|). They are also given in Table ||. 
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FIGURES 



FIG. 1. Graphical illustration of the general form of the A^-pair production process. The 
creation of a pair is described by a connected fermion line entering from and leaving to the future. 
The vacuum processes are described by all sorts of closed fermion loops. The interaction with the 
external field is shown as a cross. 

FIG. 2. One possible higher-order pair creation process, where the Pauli principle is neglected 
in the intermediate state. Interaction with the external field is shown as a cross. 

FIG. 3. Vacuum correction to the one-pair creation connected by an exchange of two lines with 
the process shown in Fig. ^. 

FIG. 4. Graphical illustration of the processes occurring in first (a) and second-order Magnus 
theory (b and c). Dotted lines denote the propagator, where only the principal part is taken, 
double lines denote electrons or positrons, which are on-shell. Interaction with the external field 
is shown as a cross. 

FIG. 5. Third order correction to the one-pair creation process in second-order Magnus theory; 
compare with Fig. ^. Double lines mean that the particles are on shell. 

FIG. 6. Lowest order (a) and the first two correction terms (b and c) to the reduced one-pair 
creation probability. 

FIG. 7. The two Feynman diagrams contributing to the one-pair creation in lowest order. The 
interaction with the external field of ion 1 and 2 is denoted by (1) and (2), respectively. 

FIG. 8. The two Feynman diagrams for the lowest-order two-pair creation process; direct term 
(a) and exchange term (b). 

FIG. 9. Comparison of the correction to the reduced one-pair creation probability with the 
lowest-order result as a function of 7. Za is set to 1. The points are the results of the calculation 
of the multiple-particle corrections, the solid line a fitted In 7 dependence. The dotted line is the 
result of the lowest-order Born approximation. 

FIG. 10. Comparison of the differential probability P{E) for 7 = 3400. Data points are the 
results of the calculation of the multiple-particle correction, the solid line results of the lowest-order 
Born calculation. Za is set to 1. 

FIG. 11. Comparison of the differential probability P{9) for 7 = 3400. Definitions are the same 
as in Fig. |l^. 
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FIG. 12. Comparison of the contributions of the direct process P(^) (dotted line) and the 
exchange process p(^) (sohd Une and data points) to the reduced two-pair creation process. Za 
is set to 1. 
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TABLES 





7 


Ion 




p{M) 




SPS 


10 


Pb 


0.63 


0.013 


2.1 % 


RHIC 


100 


Au 


1.6 


0.059 


3.7 % 


LHC 


3400 


Pb 


3.9 


0.21 


5.3 % 


LHC 


3400 


U 


6.1 


0.49 


8.1 % 



TABLE I. Comparison of the contribution of the multiple-particle correction to the re- 

duced probability with the lowest-order Born result P^^') (see Eq. (^6|)). 



n 


# Diags. 




dn 


1 


1.0000 X 10° 


1.0000000000 


1.0000000000 


3 


2.0000 X 10° 


0.3333333333 


0.6666666667 


5 


1.6000 X 10^ 


0.1333333333 


0.3777777778 


7 


2.7200 X 10^ 


0.0539682540 


0.1968253968 


9 


7.9360 X 10^ 


0.0218694885 


0.0974955908 


11 


3.5379 X 10^ 


0.0088632355 


0.0466976645 


13 


2.2368 X lO'^ 


0.0035921280 


0.0218375158 


15 


1.9038 X 10^ 


0.0014558344 


0.0100304665 



TABLE II. The number of diagrams and the coefficients Cn and d„ appearing in Magnus theory 
for different orders n. 
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